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Abstract. - It has been proved for a class of mean-field and long-range systems that the concavity 
of the thermodynamic entropy determines whether the microcanonical and canonical ensembles 
are equivalent at the level of their equilibrium states, i.e., whether they give rise to the same 
equilibrium states. Here we show that this correspondence is actually a general result of statistical 
mechanics: it holds for any many-body system for which equilibrium states can be defined and in 
principle calculated. The same correspondence applies for other dual statistical ensembles, such as 
the canonical and grand-canonical ensembles. 



Introduction. — The microcanonical ensemble (ME) 
and canonical ensemble (CE) provide two different statisti- 
cal descriptions of many-body systems at equilibrium - the 
first in terms of their energy, the second in terms of the tem- 
perature of a heat bath connected to them. Although the 
two ensembles often give compatible or equivalent predic- 
tions in the thermodynamic limit, that equivalence is not 
guaranteed by their definitions. Gibbs noticed this point as 
early as 1902 |lj and, since then, many works have shown 
that the two ensembles either give equivalent or nonequiv- 
alent thermodynamic-limit predictions depending on the 
many-body system considered or, more precisely, on the 
type of interactions acting on the particles or constituents 
of that system (see [2] for a comprehensive survey) . 

The fundamental result that has emerged from these 
works is that many-body systems involving short-range 
and stable interactions, such as short-range spin systems 
or screened Coulomb systems, always have equivalent en- 
sembles, whereas systems involving long-range interactions, 
such as gravitating systems, dipolar or unscreened Coulomb 
systems, may have nonequivalent ensembles [2j. How this 
equivalence or nonequivalence manifests itself depends on 
the level of description considered. At the thermodynamic 
level, the ME and CE become nonequivalent when the 
heat capacity is negative in the ME but is positive in the 
CE. This surprising phenomenon, which has been studied 
first in the context of gravitational systems [3jj9] , is now 
known to arise when the thermodynamic entropy of the 
ME is nonconcave as a function of the mean energy. More 
recently, it has been shown that the ME and CE can also 



be nonequivalent if they have different sets of equilibrium 
states - more precisely, if there exist equilibrium states 
in the ME, calculated at fixed energy, that do not arise 
in the CE at fixed temperature 10 . When this is the 



case, we say that the ME and CE are nonequivalent at 
the macrostate level. Various models have been shown to 
display this level of nonequivalence, including spin models 
and models of 2D turbulence 
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Given that thermodynamics is ultimately rooted in sta- 
tistical mechanics, it is natural to seek for a link between 
the thermodynamic and macrostate levels of ensemble 
equivalence, i.e., to determine whether the concavity of 
the entropy determines the macrostate equivalence of en- 
sembles. This problem was studied by Eyink and Spohn 
and later by Ellis, 
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10 , who 



Haven and Turkington 
proved that the ME and CE are, as physically expected, 
thermodynamically equivalent when they are equivalent 



at the macrostate level 17 . However, they were able to 
prove this result only for a limited class of systems that 
includes mostly mean-field-type systems. 

The aim of this paper is to show that the connection 
between the thermodynamic and macrostate levels of en- 
semble equivalence is actually a general result of statistical 
mechanics: it holds for any many-body system for which 
equilibrium states can be defined and in principle calcu- 
lated. This opens up the study of ensemble equivalence 
to many systems that do not fall in the class considered 
most critically, gravitating systems, 
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which are 

the source of much of the current research on long-range 
interactions 2,9. This general result also implies that 
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the ME and CE are always equivalent for any macrostates 
of any short-range systems having concave entropies and, 
more generally, any systems that do not have first-order 
phase transitions in the CE. 

The rigorous proofs of the results stated in this paper will 
be published elsewhere 



with 
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Here we focus on presenting 
the physical meaning of these results, as well as the main 
ideas involved in proving them, following 10 . We illustrate 



the more general nature of our results by discussing the 
equivalence of the ME and CE for a simple spin model not 
completely covered by the results of 10 . We also derive, 



from our equivalence results, general results about phase 
coexistence and first-order phase transitions. Finally, we 
discuss their generalization to ensembles other than the 
ME and CE. 

Framework. — We consider a general classical N- 
particle system described by the Hamiltonian or energy 
function H(oj), where ui € A denotes the microstate of that 
system and A the microstate space. The mean energy or 
energy per particle is denoted by h(u) = H(u)/N. We 
also consider a macrostate of that system, defined math- 
ematically as a function M(oS) of the microstates. This 
macrostate can represent, for example, the magnetization 
of a spin system, the empirical distribution of velocities or 
positions of a gas of N particles, or a combination of any 
observables defined similarly as functions of uj. 

The problem that we are concerned with is to define 
the set of equilibrium values of M in both the ME and 
CE in the thermodynamic limit, and to determine whether 
these sets are equivalent, i.e., whether they can be put in 
some correspondence. Mathematically, the microcanonical 
equilibrium values of M, which we denote collectively as 
£ are identified as the typical values (viz., global maxima 
or concentration points) of the ME probability distribution 
P U {M = m) obtained by conditioning a uniform (prior) 
probability distribution P on A on the "constrained" set 
of microstates such that h(u>) = u. In symbols, this is 
expressed as 



P u (m) = P(M = m\h = u) 



P(M = m,h = u) 
P{h = u) ' 



(1) 



Note that, for simplicity, we use the same P to denote the 
prior distribution P(ui) on A and probabilities in general. 
Thus P(M = m) denotes the probability (density) of the 
event M = m calculated with respect to P(w). 

In a similar way, the canonical equilibrium values of M, 
denoted by £p, are identified as the typical values of the CE 
probability distribution Pp{M = m) of M obtained by in- 
tegrating the Gibbs distribution over the set of microstates 
such that M{ui) = m. Thus, 



Pp{M = m) = [ Pp(u)6u, 



where 



(2) 



(3) 



Z{fi) 



(4) 



as the partition function. The parameterization of the 
different distributions reflects of course the parameters 
defining each ensemble: the mean energy u (or the energy 
Nu) for the ME, and the inverse temperature j3 (or tem- 
perature T = {ksP)~ l with fee the Boltzmann constant) 
for the CE. 

Mixture of ensembles. — It is obvious from the con- 
struction and physical interpretation of the ME and CE 
that P u (m) and Pp(m) are in general two very different 
probability distributions of the same macrostate. The point 
about equivalence of ensembles, however, is that the typi- 
cal values of M determined by these distributions, which 
define respectively £ u and £p, may be equivalent in the 
sense thatm for a given u, there might be a given f3 such 
that £ u = £p. The essential idea involved in establishing 
this equivalence, which can be traced back to Gibbs [l], is 
that the CE is a probabilistic mixture of MEs. 

To explain the meaning of this last statement, simply 
consider the Gibbs canonical distribution of Eq. ([3| . Since 
this distribution depends only on (3 and H(oj), it is clear 
that all microstates u> having the same energy have the 
same probabilistic weight, which implies that the condi- 
tional probability distribution Pp(us\h = u), obtained by 
conditioning Pp{uj) on the set of microstates oj such that 
h(cj) = u, must be uniform over that constrained set of 
microstates. Since P u (u>) has exactly this property, we 
must therefore have Pr(uj\u) — P u (lu) for all to £ A and 
by extension Pp(m \u) = P w (m). 

With this result, we can use Bayes's Theorem to write 

Pp{m)= / P p (m\u) P p (u) du = / P u {m)P f) {u)du, 



(5) 

where Pp(u) is the probability distribution of the mean 
energy h in the CE. In terms of the density of states fi(u), 
this distribution is simply given by 
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Z(fi) 



(6) 



The second equality obtained in ([5| shows that the CE is 
a superposition of MEs weighted by the canonical prob- 
ability distribution Pp (u) of the mean energy. It is this 
superposition, which is valid for all N, that we refer to as 
a probabilistic mixture of MEs. 

Results. — The idea that the CE is a superposition of 
MEs is simple but far reaching because it implies, by taking 
the thermodynamic limit, that the concentration points of 
Pp{m) are the concentration points of P u (m) where Pp{u) 
concentrates. In other words, £p must be given by all the 
sets £" such that u is an equilibrium value of the mean 
energy in the CE at inverse temperature /3. To put this 
statement in mathematical form, let us denote by Up the 
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set of these mean energies realized at equilibrium in the 
CE at (3. Then we obtain from Eq. 



(7) 



This set covering result is a central result of this paper. 
The equivalence or nonequivalence of the ME and CE 
follows from this result by determining whether Up contains 
one or more equilibrium mean energies u and whether a 
given u is a member of Up for some /3 S K. These questions 
are known to be answered [17] by the concavity properties 
of the microcanonical thermodynamic entropy or entropy 
density, defined by the usual limit 



s(u) = lim — lnJ7(u). 



(8) 



We state next the essence of the relationship between the 
concavity of s(u) and the form of Up, and then state the 
consequence of this relationship for the equivalence of the 
ME and CE as derived from Eq. ([5]). For the definition of 
concave points of s(u), see [17]. A basic reference guide to 
this definition, which is sufficient for the purpose of this 
paper, is given in Fig. [TJ 

As in 1 10 , we are led to consider three cases: 
Case 1 (Equivalence): If the entropy s is strictly 
concave at u (point a in Fig. [TJ), then there exists an 
inverse temperature /? € M such that Up is a singleton set 
containing only u. For this inverse temperature, Eq. ^ 
then reduces to £ u = £p. In this case, the equilibrium 
states of M in the ME are all realized as equilibrium states 
in the CE. 

Case 2 (Nonequivalence): If s is nonconcave at u 
(point b in Fig. [TJ, then u ^ Up for all /3 g R, so that 
£ u ^ £p for all j3 € R. In fact, it can be proved in this 
case that £ u f\£p = % for all /3 e R. 

The physical interpretation of this case is that the equi- 
librium states obtained in the ME for values of u where s(u) 
is nonconcave are not realized as equilibrium states in the 



CE for any temperature 17 . We call these states micro- 



canonical nonequivalent states. Models having such states 
include the mean-field Potts model 13 , the mean-field 
Blume-Emery-Griffiths model [12], and the point- vertex 
model of 2D turbulence [14] . 

Case 3 (Partial equivalence): If s is concave but not 
strictly concave at u (point c in Fig. [TJ), then u <EUp for 
some /3 € R, but u is not the only member of Up, which 
implies from Eq. ^ that £ u C £p. 

This last result follows because £p in this case is made up 
of different microcanonical equilibrium states at different 
mean energies, one of which is u. Thus £p contains all the 
equilibrium states of the ME at u but may also contain 
more equilibrium states not seen in £ w ; hence the name 
"partial equivalence" . If any of the many sets that compose 
£p is different from £ u , which is expected since different 
energies should give rise to different equilibrium states, 
then the previous result is strengthened to £ u C £p. That 
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Fig. 1: Microcanonical entropy function s(u) (black) and its 
concave envelope s**(u) (grey), (a) Strictly concave points of 
s; (b) nonconcave points of s; (c) Concave points of s that are 
not strictly concave. 



is, in this case £p does contain equilibrium states that are 
not contained in £ u . 

Relationship with phase transitions. Partial 
equivalence of ensembles is a boundary case between the 
equivalent and nonequivalent cases, which corresponds 
physically to a phase coexistence accompanying a first- 
order phase transition in the CE [17]. To understand this 
point, consider the two points labelled by c in Fig. [T] These 
two points are not strictly concave because they share the 
same supporting line, which corresponds to the line of the 
concave envelope s**(u) of s(u). Let us denote the mean 
energies of these points by ui and Uh, respectively, and let 
/3 C be the slope of their supporting line. The non-strict 
concavity of these points then implies that Pp a (u) is a 
bimodal distribution with two equal-height peaks at ui and 
~, so that Up a = {ui,u h } and £p c = £ u ' U£ Uh . This 
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Uh _ 

shows that the CE at /3 C is composed of many equilibrium 
states of the ME. If we regard each of these microcanonical 
sets as a "pure phase" of the ME, then we can say that 
the CE shows phase coexistence at j3 c . 

This result is consistent with the fact that the CE shows 
a first-order phase transition at j3 c - Indeed, it is known 
that if the entropy s(u) is nonconcave, for example as in 
Fig. [lj then the canonical free energy, 



¥>(/?)= lim -^lnZ(p), 



(9) 



is nondifferentiable at some point, corresponding in Fig. [TJ 
to (3 C . We thus see that the set of equilibrium states of 
any system in the CE is, at a point of first-order phase 
transition in the CE, the combination of two or more sets 
or "phases" of equilibrium states of the ME. This does 
not imply that each equilibrium state in the CE is itself 
a superposition or a mixture of two or more equilibrium 
states of the ME; only that the whole set of equilibrium 
states in the CE is a such a mixture. 

Note, finally, that the connection between nonconcave 
entropies and nondifferentiable free energies also goes in 
the opposite direction, in that, if ip(/3) is everywhere dif- 
ferentiable, then s(u) is necessarily strictly concave. We 
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therefore obtain as a general physical result that the ME 
and CE are equivalent at the macrostate level for any 
systems that do not have first-order phase transitions. 

Comparison with previous results. The result of 
Eq. ([5]) and the three cases stated above were first obtained, 



as mentioned before, by Ellis, Haven and Turkington f 



The contribution of the present paper is to re-derive these 
results under much more general conditions than originally 
done. All that is required here is for the distributions 
P u (m) and -Pg(m) to exist, so that the equilibrium sets £ u 
and £r can be constructed. This is a very weak assumption, 



which amounts essentially to saying that the results of 10 



actually hold provided that the thermodynamic limit of 
the ME and CE is well defined. In fact, it can be proved 
that the existence of £ u and £p is essentially equivalent to 
the existence of s{u) and <p(/3), respectively, so that our 
results about the equivalence and noncquivalence of the 
ME and CE simply hold provided these thermodynamic 
functions exist. 

In 10 , the same results were obtained by comparison 
under two technical conditions. The first is that the mean 
energy function h(uj) can be rewritten as a function of some 
macrostate M{uS) in the thermodynamic limit. Mathemat- 
ically, this amounts to saying that there exists a function 
h(m) of the macrostate M, called the energy representation 
function, such that 



lim \h(u) 

iV— >-00 



h(M(oj))\ = 



(10) 



uniformly for all uj G A. The second condition is that there 
exists an entropy function s{m) associated with M. This 
macrostate entropy is defined similarly as s{u) by replacing 
the density of states fi(u) with the density of states fi(m) 
associated with the microstates having a given macrostate 
value M(ui) — m (see [lOj for the precise definition). 

These two conditions limit greatly the type of systems for 
which the equivalence of the ME and CE could be analyzed. 
The first condition in particular is akin to a mean-field 
requirement, known to be satisfied for mean-field systems 
and certain long-range interacting systems 10 12 19 . It 



does not apply in a uniform sense to gravitational systems, 
among other long-range systems, and certainly not for 
systems with short-range interactions. Moreover, even for 
mean-field or long-range systems, there is typically only one 
macrostate that qualifies as a macrostate for which h and 
s exist, which means that the equivalence of the ME and 
CE can only be discussed for this particular macrostate. 

In our version of the equivalence and nonequivalence 
results, these restrictions are completely lifted: we can 
discuss the equivalence or nonequivalence of the ME and 
CE for any systems and any macrostates or observables 
provided that the equilibrium values of these observables 
can be defined in each ensemble in the thermodynamic 
limit. When this is the case, the ME and CE are equivalent 
essentially if and only if the microcanonical entropy is 
concave in the thermodynamic limit as a function of the 
mean energy. 



Application. — To illustrate the generality of our re- 
sults, compared to those of u0\, we now discuss the equiv- 
alence of the ME and CE for the a-Ising model, defined 
by the Hamiltonian 



H = 



J 



N 1 



N 

E 

i>j=i 



1 SiSj 

" -j\ a 



(ii) 



where Si = ±1 is a spin variable. This model has been 
extensively studied because of its rich phase transition 



behavior 20|-|23 , which includes a suspected Kosterlitz- 

For J > and 



Thouless transition for a 
< a < 1, Barre et 



2 24 25 



19 were able to show that the 



model admits an energy representation function h and 
macrostate entropy s for a particular macrostate, namely, 



the macroscopic magnetization profile m{x) (see 19 for 



the precise definition of this macrostate). Given that the 
entropy s(u) is concave as a function of u for these pa- 
rameters, they then concluded that the ME and CE are 
equivalent at the level of m{x). 

For a > 1, there is no known energy representation 
function, so the equivalence of the ME and CE cannot be 
discussed for this case using the results of 10 . However, 



the entropy s(u) exists for these parameters and is every- 
where concave, since it is known that the model does not 
show any first-order phase transitions. From our results, 
we therefore conclude that the ME and CE are equivalent 
for this model at the level of m(x). More importantly, they 
are also equivalent at the level of any macrostate, such as 
the standard magnetization or the spin distribution, which 
do not satisfy the requirements of 10 . Similar conclusions 



can be reached for other models and macrostates, includ- 
ing short-range models, such as the Ising model at the 
level of its magnetization, and long-range models such as 
self-gravitating particle models [18] . 

Idea of the proof. — We conclude by giving an 
overview of how the results of this paper are proved at a 
rigorous level. The first step, following [To], is to use large 
deviation theory 26 to derive a large deviation principle 
(LDP) for each of the distributions P u (m) and Pp(m), and 
to use these LDPs to define (in the thermodynamic limit) 
the sets £ u and £p of equilibrium values of the macrostate 
M in the ME and CE, respectively. 

Heuristically, we express these LDPs by the approxima- 
tions 

-NI u (m) 



and 



P u {m) 



Pp{m) 



-NIf,{m) 



(12) 



(13) 



which are assumed to hold with sub-exponential correction 
terms in the limit N — > 00. Given the rate functions I u (m) 
and Ip(m), the sets £ u and £p are defined, respectively, as 
the sets of global minima of I u (m) and Ip{m) 
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Next, we use the idea that the CE is a probabilistic 
mixture of MEs, i.e., Eq. to obtain the following large 
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deviation approximation: 



Pp{m) 



= -JV{7"(m) + ./ f) (u)} du 



which involves the rate function I u (m) of P u (m) and the 
rate function Jp (it) associated with the LDP of the distri- 
bution Pp{u) of the mean energy h in the CE: 



(14) I wish to thank R. Nerattini, C. Nardini, G. L. Sewell 
and R. J. Harris for useful comments. This work was 
partially supported by an Interdisciplinary Fellowship from 
RCUK. 
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